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Abstract. We investigate the nonlinearity-assisted quantum tunneling and formation 
of nonlinear collective excitations in a matter-wave interferometer, which is realised by 
the adiabatic transformation of a double-well potential into a single-well harmonic 
trap. In contrast to the linear quantum tunneling induced by the crossing (or avoided 
crossing) of neighbouring energy levels, the quantum tunneling between different 
nonlinear eigenstates is assisted by the nonlinear mean-field interaction. When the 
barrier between the wells decreases, the mean-field interaction aids quantum tunneling 
between the ground and excited nonlinear eigenstates. The resulting non- adiabatic 
evolution depends on the input states. The tunneling process leads to the generation 
of dark solitons, and the number of the generated dark solitons is highly sensitive to 
the matter-wave nonlinearity. The results of the numerical simulations of the matter- 
wave dynamics are successfully interpreted with a coupled-mode theory for multiple 
nonlinear eigenstates. 
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1. Introduction 

Matter-wave interferometry involves coherent manipulation of the external or internal 
degrees of freedom of massive particles lU [2] . Utilizing the well-developed techniques of 
trapping and cooling, the matter-wave interferometers have been realised with atomic 
Bose- Einstein condensates (BECs) [31 SI [5], E] . Almost all BEC interferometers based on 
spatial interference measure the phase coherence by merging two initially separated 
condensates [31 HI El E]. To recombine two condensates confined in a double- well 
potential, one has to transform the double-well potential into a single-well potential 
by decreasing the barrier height. 

Intrinsic interparticle interactions in atomic condensates have stimulated various 
studies of the nonlinear behaviour of condensed atoms [7]. A balance between matter- 
wave dispersion and nonlinear interaction supports a number of nontrivial collective 
excitations, including bright solitons in condensates with attractive interactions [H [9] 
and dark solitons in condensates with repulsive interactions [IHl |TT] . In a harmonically 
trapped condensate with repulsive interparticle interactions, the nodes of excited 
nonlinear eigenstates correspond to dark solitons [16], so that the formation of dark 
solitons can be associated with populating excited states [El [17]. Several methods 
of condensate excitation have been suggested, the most experimentally appealing ones 
involving time-dependent modifications of trapping potentials [18]. The operation of 
BEC interferometers and splitters based on spatiotemporal Y- and X-junctions [121 E] 
is greatly affected by the possibility of nonlinear excitations. The nonlinear excitations 
in BEC interferometers with repulsive interparticle interactions lead to the generation of 
dark solitons [131 [H], and can be utilised to enhance the phase sensitivity of the devices 
[HiTS]. 

The extensively explored mechanisms for population transfer between different 
eigenstates of a trapped BEC include non-adiabatic processes [12], Josephson tunneling 
[T9l [20] . and Landau- Zener tunnehng [211 l22l [23], which are also responsible for 
population transfer in linear systems. However, in a sharp contrast to linear systems, 
the quantum tunneling between different nonlinear eigenstates can be assisted by the 
nonlinear mean-field interaction even in the absence of crossing (and avoided crossing) 
of the energy levels. Up to now, this peculiar type of quantum tunneling remains poorly 
explored. 

In this paper, we explore the intrinsic mechanism for the quantum tunneling 
assisted by repulsive nonlinear mean-field interactions in a matter-wave interferometer. 
We consider the dynamical recombination process of a BEC interferometer, in which 
an initially deep one-dimensional (ID) double- well potential is slowly transformed 
into a single-well harmonic trap. Our numerical simulations, employing a time- 
dependent ID mean-field Gross-Pitaevskii (CP) equation, show that multiple moving 
dark solitons are generated as a result of the nonlinearity-assisted quantum tunneling 
between the ground and excited nonlinear eigenstates of the system, and the qualitative 
mechanism is independent on the particular shape of the symmetric double-well 
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potential. Furthermore, the number of the generated dark sohtons is found to be 
highly sensitive to the strength of the effective nonlinearity that in turn depends on 
the total number of condensed atoms and the atom-atom s-wave scattering length. The 
population transfer between different nonlinear eigenstates caused by the nonlinearity- 
assisted quantum tunneling can be quantified by a coupled-mode theory for multiple 
nonlinear eigenstates of the system. 




time 



Figure 1. Schematic diagram. Top: time-dependence of the barrier height B{t). 
Bottom left: initial BEC density distribution (shaded) at t = in the double- well 
potential (solid line). Bottom right: density distribution (shaded) at i = 80 in the 
single- well potential (solid line). 



2. Model and numerical results 

We consider a condensate under strong transverse confinement, muj W + z'^)/2, so that 
the 3D mean-field model can be reduced to the following ID model \24\ : 
d 

zh—^{x,t) = Ho^{x,t) + X\^{x,t)f^{x,t), (1) 

where Hq = -{h^ /2m){d'^ / dx^) + V{x,t), m is the atomic mass, A > characterizes 
the effective nonlinearity which we assume to be repulsive, and V{x,t) is an external 
potential. If the condensate order parameter \l/(x,t) is normalised to one, the effective 
nonlinearity A = 2NasUjph is determined by the total number of atoms N , the s-wave 
scattering length and the transverse trapping frequency ujp [25]. In what follows we 
use the dimensionless version of the model equation obtained by choosing the natural 
units of m = ^ = 1. 

We assume the time-dependent potential V{x, t) as a spatiotemporal Y-shape 
potential generated by the superposition of a ID time-independent harmonic potential 
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Figure 2. Evolution of the condensate density (left) and phase (right) of the ground 
and first-excited states for the system of the effective nonlinearity A = 20. The first 
and second rows correspond to the evolution of the ground and first-excited states, 
respectively. 



and a time-dependent Gaussian barrier (see Fig. 1): 
V{x, t) = V + B{t) ■ exp [ J , 



(2) 



where u is the trapping frequency, d is the barrier width, and the barrier height depends 
on time as follows: 

fBo-at, for t<Bo/a, 
^^'^ = \0, for t>Bo/a, 

where a is the rate at which the barrier between the wells is ramped down. When the 
barrier height B{t) > cu^rf^, the time-dependent potential is a double- well potential with 
two minima at a; = ±d^j2 ln[i?o/ (rf^tu^)]. Thus, the ID description is valid for weak 

longitudinal confinements satisfying \n[Bo/{d'^uj'^)] <ti ujp and u <^ ujp. To 

ensure the adiabatic evolutions of the symmetric and antisymmetric initial eigenstates, 
the rate a must be sufficiently small. In Fig. 2, we show the evolutions of the ground and 
first-excited eigenstates of the system of the effective nonlinearity A = 20 and ramping 
rate a = 1/4. For such a small value of a, both the ground and first excited eigenstates 
of the initial double-well potential adiabatically evolve into the corresponding ground 
and first excited eigenstates of the final single-well potential. This means that the 
non- adiabatic effects are negligible for such a slowly varying process. 
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The usual double-well BEC interferometers involve condensates trapped in a 
double-well potential before recombination. Below, we consider the case of the initial 
state of the BEC being fully localised in a single well of a symmetric double-well potential 
with a sufficiently high barrier, so that there is no significant overlap between the 
Wannier states of the two wells. The fully localised initial state can be viewed as 
the equal-probability superposition of the ground and first-excited eigenstates, so that 
it can be used to observe the interference of these two eigenstates [26]. For BECs 
trapped in such a deep potential, the mean-field ground and first excited states are 
degenerate or quasi-degenerate. Even for low barriers, if the tight-binding condition 
is still satisfied, the two-mode approximation will give the picture of a classical Bose- 
Josephson junction. In the framework of second quantisation, the system obeys a two- 
site Bose-Hubbard Hamiltonian. In this fully quantum picture, a completely localised 
initial state corresponds to the highest excited state for repulsive interactions. This 
state exhibits degeneracy with the sub-highest excited state [27], which corresponds to 
the bistability in a classical Bose-Josephson junction [28] . 
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Figure 3. Evolution of the condensate density (left) and phase (right) for different 
values of the effective nonlinearity A. Cases (a-d) correspond to A = 0, 5, 10 and 15, 
respectively. 
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Since initially there is no overlap between the two Wannier states, the quantum 
tunneling between those states is negligible. As the barrier height gradually decreases, 
the overlap between two Wannier states becomes more significant. Then both the 
quasi- degeneracy between the ground and first excited states in the mean-field picture 
and the quasi- degeneracy between the highest excited and sub-highest excited states in 
the quantum picture break down. The quantum tunneling of the fully localised state 
becomes more pronounced as the barrier height is decreasing. Due to the very slow 
reduction in the barrier height, the kinetic energy stays small during the whole process, 
and the quantum tunnelling dominates the dynamics. However, the over-barrier hopping 
could occur in a rapidly varying process, in which case the kinetic energy can exceed 
the potential barrier. 

To explore the dynamic evolution, we numerically integrate the GP equation with 
the well-developed operator-splitting procedure and the absorbing boundary conditions. 
In Fig. 3 we show the time evolution of the condensate density and phase for the trapping 
frequency cu = 0.27r, the initial barrier height Bq = 20.0, the barrier width d = the 
ramping rate a = 1/4, and different values of the effective nonlinearity A. For the chosen 
small ramping down rate, all symmetric (antisymmetric) states of the deep double-well 
potential will adiabatically evolve into the corresponding ground (or excited) states of 
the single-well potential. 




Figure 4. Formation of dark solitons in tlie system witli the effective nonlinearity 
A = 15. Left: density distributions |^'(a;,i)p for different times. Right: phase 
distributions (fi^x, t) for the corresponding density distributions. 
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Evolution of the fully localised initial state strongly depends on the values of the 
effective nonlinearity A. In the linear case (A = 0), the fully localised initial state can 
be viewed as an equal-probability superposition of the symmetric and antisymmetric 
states, so that, according to the adiabatic theorem, the evolving state is always the 
equal-probability superposition of the ground and the first excited states of the system. 
Due to the nonlinear interactions, the superposition principle becomes invalid, and 
the resulting behavior can be interpreted as the coupled dynamics of the ground and 
multiple excited states of the nonlinear system. Akin to the linear systems, the quantum 
tunneling appears once the quasi-degeneracy between the ground and the first excited 
state is broken, and gradually becomes significant with decreasing barrier height. The 
time scale on which the quantum tunneling appears in the nonlinear system shortens 
with the growth of nonlinear interaction strength A. The excited nonlinear states of the 
BEC in a single-well potential can be thought as stationary configurations of single or 
multiple dark solitons [16] . As a result of the population transfer to such excited modes, 
the condensate develops multi-peak distribution with significant phase gradients across 
density notches between neighboring peaks. These notches are dark or gray solitons 
with well-defined phase gradients close to vr (see Fig. 4). 

The number of dark solitons formed in this process varies with the effective 
nonlinearity A (see Fig. 5). This dependence exhibits multiple plateaus as the effective 
nonlinearity A changes, as shown in Fig. 5. Given the relationship between the 
nonlinear interaction strength and the key parameters of the system, A = 2NasOOph, 
one can control the number of generated solitons by adjusting the s-wave scattering 
length with the Feshbach resonance, the total number of atoms in the condensate with 
initial preparation, and/or the transverse trapping frequency by tuning the transverse 
trapping field strength. The number of solitons remains unchanged for a long period 
of time before multiple collisions between solitons take place. The inelastic collisions 
lead to the radiation of small-amplitude waves, and after a large number of collisions 
the number of solitons oscillating in the trap changes. 
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Figure 5. Number of dark solitons generated in the condensate at i = 80 versus the 
effective nonlinearity parameter A. 
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3. Modal decomposition 

To obtain the quantitative picture of the population transfer, we decompose an arbitrary 
state of our time-dependent system as [15^ ?], 

TV 

vI/(x,t)=5:C,(t)0,(a;,t), (4) 

j 

where (j)j{x,t) is the j-th stationary state for the nonlinear system with the potential 
V{x,t), which obeys the (dimensionless) equation: 



lXj{t)(l)j{x,t) 



2 



^j{x,t) + \(j)%x,t). (5) 



Here ^j{t) is the chemical potential for the j-th stationary state. For each instantaneous 
form of the time-dependent potential, the nonlinear eigenstates {(j)j{x,t)} form an 
orthogonal set, similarly to their linear counterparts [151 [16], [29] . Due to the nonlinear 
interparticle interactions, there exist additional stationary states (e.g., self-trapped 
states) which have no linear counterparts [291 ISO]- Nevertheless, every stationary 
state of the nonlinear system can be composed from the orthogonal basis {(f)j{x,t)}. 
The nonlinear eigenstates of time-independent potentials are also time-independent. 
However, due to the violation of the superposition principle, population transfer between 
different nonlinear eigenstates also occurs in time- independent systems |29j. This type 
of population transfer, which originates from the exchange collisions between atoms in 
different eigenstates, will be considered in detail below. 

The population dynamics for different nonlinear eigenstates can be described by 
the evolution of the complex coefficients Cj{t), which obey a series of coupled first-order 
differential equations. 



JCi{t) ^ ^ 



dt 



E'o' +j:Qk,k'Cmck'it) 



C,{t). (6) 



k,k' 

Due to the conservation of the total number of particles, Cj{t) satisfy the normalisation 
condition J2j 1^(^)1^ = 1- Here the linear coupling parameters are 

= / 0r(x,t)ffo0,(x,t)rfa;, (7) 

and the nonlinear coupling parameters are 

Qfeffe'W = ^ / (pKx^'t)<P*kix,t)(j)kix,t)(j)j{x,t)dx. (8) 



For a spatially symmetric potential V{x,t) = V{—x,t), we have Q\^^ki{t) = 0, and then 
{k + k' + I + j) are odd integer numbers. 

In our numerical simulations, we generalise the direct relaxation method for linear 
quantum systems [31j to calculate the eigenstates and their eigenvalues (chemical 
potentials) for our nonlinear system with different effective nonlinearities at any moment 
of time. Projecting the condensate wavefunction \&(a;, t) onto the nonlinear eigenstates 
(l)j{x,t), we find the population probabilities Pj(t) = |Cj(t)p = |/ 0*(x, t)\E'(x, t)(ixp 
which depend on the time and the effective nonlinearity A. 
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Figure 6. (a) Time evolution of population probabilities in different eigenstates for 
A = 2, where Pq and Pj {j — 1,2,3) are the population probabilities of the ground 
state and the j— th excited state, respectively. The total probability of the first four 
eigenstates is denoted as Ptot = Po + Pi + P2 + Pa- (b) The corresponding population 
evolution for i > 80 obtained from the coupled-mode equation (6) with first four lowest 
eigenstates. (c) Population probabilities for the first six lowest nonlinear eigenstates 
at t = 80. 



In Fig. 6(a), we show the time evolution of the population probabilities Pj{t) for 
the effective nonlinearity A = 2.0. Here we only consider four lowest eigenstates (i.e., 
iV = 4), so that Pq is the ground state population probability, Pj (j = 1,2,3) are the 
population probabilities of the j— th excited state, and Ptot = P0+P1+P2+P3 is the total 
probability of the first four eigenstates. For t < 50, the population probabilities keep 
almost unchanged. In the region of 50 < t < 80, we observe a fast population transfer 
from the ground state to the first excited state. After the recombination of the two wells, 
i.e. for t > 80, the populations in different nonlinear eigenstates oscillate with time, 
even though the nonlinear system has a time-independent potential, time-independent 
eigenstates, and no energy degeneracy between neighbouring eigenstates. This behavior 
differs drastically from the linear dynamics where populations in different eigenstates 
always remain unchanged. We find that at this value of the effective nonlinearity the 
total population probability Ptot(^) in the first four eigenstates is always close to one. 
The low-frequency population oscillations are dominated by the linear coupling between 
different modes and the high-frequency ones are due to the nonlinear cross-coupling of 
the nonlinear modes which corresponds to the exchange collision of atoms in different 
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eigenstates. 

For small A, the dynamics of Pj after the merging of the two wells can be 
approximately captured by the projection of the BEC state at the moment of the 
merging (here t = 80) onto the set of stationary nonlinear states (pjlx) of the single- 
well potential Vo{x) of B{t) = 0. This is confirmed in Fig. 6(b), where we employ the 
coupled-mode theory ([6]) with = 4 eigenstates of Vo{x) [cf. Fig. 6(a)]. The number of 
eigenstates, N, that must be considered in the coupled-mode theory, increases with the 
effective nonlinearity. The highest-order mode (j = A^) of the harmonic potential Vo{x) 
with significant (non-zero) excitation probability P/v at the merging time will therefore 
determine the number of dark solitons that are likely to be formed. Figure 6(c) shows 
the excitation probabilities at t = 80 for the first six lowest eigenstates of Vq{x) and 
different A. By comparing the number of significantly excited states for different values 
of A, one can see that the number of dark solitons formed is indeed approximately 
determined by the highest-excited nonlinear mode of the harmonic trap that is still 
sufficiently populated. For instance, A^ = 1 solitons are expected to form for A = 2, and 
AT = 2 for A = 10 (cf. Fig. 5). 

4. Conclusions 

We have explored the nonlinearity-assisted quantum tunneling and formation of 
nonlinear collective excitations in the matter-wave interferometer based on a time- 
dependent double-well potential, dynamically reconfigured to form a single-well 
harmonic trap . In contrast to the Josephson tunneling and Landau-Zener tunneling, the 
nonlinearity-assisted quantum tunneling is brought about by the nonlinear inter-mode 
population exchange scattering. The excitations caused by this type of tunneling lead 
to the dark soliton generation in the process that differs dramatically from the phase 
imprinting [101 [11] or condensates collisions [13]. The number of generated solitons can 
serve as a sensitive measure of the degree of the nonlinearity in the system. With the 
well-developed techniques for preparing and manipulating condensed atoms in double- 
well potentials [31 HI [5], [6] , loading the condensed atoms in one well of a deep double- well 
potential and adjusting the barrier height, the experimental observation of this effect 
seems feasible. 

Note added: After this manuscript was prepared for submission, the group of Peter 
Engels from the Washington State University reported the experimental observation of 
matter- wave dark solitons due to quantum tunneling [32], in the process of sweeping a 
high potential barrier from one edge of the trap to the other. 

The authors thank R. Gati and M. Oberthaler for stimulating discussions. This 
work was supported by the Australian Research Council (ARC). 
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